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Abstract 

The current paper is a corrected version of our previous paper (Olami et al., PRE 
55 (3), (1997)). Similarly to previous version we investigate the problem of flame 
propagation. This problem is studied as an example of unstable fronts that wrinkle 
on many scales. The analytic tool of pole expansion in the complex plane is em- 
ployed to address the interaction of the unstable growth process with random initial 
conditions and perturbations. We argue that the effect of random noise is immense 
and that it can never be neglected in sufficiently large systems. We present simu- 
lations that lead to scaling laws for the velocity and acceleration of the front as a 
function of the system size and the level of noise, and analytic arguments that ex- 
plain these results in terms of the noisy pole dynamics. This version corrects some 
very critical errors made in (Olami et al., PRE 55 (3), (1997)) and makes more 
detailed description of excess number of poles in system , number of poles that 
appear in the system in unit of time, life time of pole. It allows us to understand 
more correctly dependence of the system parameters on noise than in (Olami et al., 
PRE 55 (3),(1997)) 

Keywords: flame propagation, pole dynamics, unstable front, random noise, self- 
acceleration 

1 Introduction 

This article considers the very interesting problem of describing the nonlinear stage of 
development of hydrodynamic instability of the flame. This problem can be considered 
for ID (channel propagation) , 2D (cylindrical case) , 3D (spherical case). The direct 
numerical simulations can be made based on of the NavierStokes equations including 
chemical kinetics in the form of the Arrhenius law 0]. For 2D and 3D cases we can 
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see experimentally observable effects 0, of the self-acceleration of the front of a 
divergent flame, the formation of cellular structure, and other effects. 

Much more simple equation for flame front propagation can be obtained 0,0. 
It is Michelson-Sivashinsky approximation model. The Michelson-Sivashinsky model 
assumes very serious limitations, such as the smallness of the coefficient of gas ex- 
pansion and, consequently, the potential flow in the combustion products and fresh 
mixture, weak nonlinearity, the assumptions of the stabilizing effect of the curvature 
of the flame front and a linear dependence on the curvature of the normal speed, and 
others. The calculations with a noise term was performed in both one-dimensional and 
two-dimensional formulations of the problem in our previous papers 0, 0, El, (2l> 
1 10 1, 1 1 1 1 and recently by Karlin and Sivashinsky fl2l , lfl3l for ID, 2D and 3D cases. 

Interest of this model stems, firstly, from the fact that despite the serious limita- 
tions, this model can qualitatively describe the scenario of hydrodynamic instability 
and, in particular, the self-acceleration of the front of a divergent flame, the formation 
of a cellular structure, and other effects. Secondary, this nonlinear model has exact 
solutions which can be constructed on the basis of pole expansions [14], fl5l . fl6l . 
ifTTl . This pole expansion method got development in the following papers. There was 
investigated relationships between pole solutions and partial decomposition in Fourier 
series |T8ll , between the gas flow field and pole solutions [ 19 1. The future development 
can be found in our previous papers and the correspondent references inside of |6|, |7|, 

@, 0, ED, HE 

It must be mentioned that the simplest ID case of the flame front propagation is 
very important. It was the main reason for creation this papers considering in detail 
the ID case. Such investigation of this case allows us to understand qualitatively and 
quantitatively the pole dynamics. This understanding is a basis for the consideration 
more complex 2D and 3D cases. The our paper |8| clearly demonstrate this fact. The 
cellular structure, acceleration exponent for 2D case was found on the basis of ID 
results. 

The future development of our results is gotten in papers [20|, |2TI . Il22l . In the 
very interesting paper |20| the noise term was considered in the poles-like form. The 
numerical results for ID case are in a good consistence with the theoretical results of 
this paper. 

This paper is update version of our previous version |6 1. In this version we correct 
some error. For example, we choose incorrect spectrum width of white noise with con- 
stant amplitude during changing size of channel L. As result we didn't obtain in |6| 
saturation of the frame front velocity during increasing L in contradiction with this pa- 
per. The second important error connected to number of regimes defined on the "phase 
diagram" of the system as a function of L and /. Indeed not three but four such regimes 
exist. The noise influence theory developed in [6 1 described transition from regime I 
to regime II and describe only these two regimes. The regime II can not be observed 
because of numerical noise. The observable regime III can not be explained by the help 
theory developed in [6 |. To make such explanation we calculate numerically and try 
explain analytically such values as excess number of poles in system , number of poles 
that appear in the system in unit of time, life time of pole ATI . The developed theory 
can explain the existence of small dependence parameters of problem on the noise in 
regime III. 
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The rest of the paper is organized as follows. We begin by presenting equations 
of motion and pole-decomposition in the channel geometry(Section[2]l.Next Section[3] 
describes acceleration of the flame front, pole dynamics and noise. And finally (Sec- 
tion |4|i we give summary and conclusions. 



2 Equations of Motion and Pole-decomposition in the 
Channel Geometry 

It is known that planar flames freely propagating through initially motionless homoge- 
neous combustible mixtures are intrinsically unstable. It was reported that such flames 
develop characteristic structures which include cusps, and that under usual experimen- 
tal conditions the flame front accelerates as time goes on. A model in 1 + 1 dimensions 
that pertains to the propagation of flame fronts in channels of width L was proposed in 
0. It is written in terms of position h(x, t) of the flame front above the x-axis. After 
appropriate rescalings it takes the form: 



dh(x,t) _ 1 
dt ~~ 2 



dh(x, t) 



d 2 h(x,t) 
dx 2 



+ I{h(x,t)} + 1 



(1) 



The domain is < x < L, v is a parameter and we use periodic boundary conditions. 
The functional I[h(x,t)] is the Hilbert transform of derivative which is conveniently 
defined in terms of the spatial Fourier transform 



h(x,t) = I e ikx h{k,t)dk 



I[h(k,t)] = \k\h(k,t) 



(2) 
(3) 



For the purpose of introducing the pole-decomposition it is convenient to rescale the 
domain to < 9 < 2ir. Performing this rescaling and denoting the resulting quantities 
with the same notation we have 



dh(6,t) 
dt 



1 

2L 2 



dh(0,t) 



~i{h(e,t)} 



i 



v d 2 h(9,t) 
L 2 dO 2 



(4) 



In this equation L = L/2ir. Next we change variables to u(0,t) = dh(0,t)/d6. We 
find 

du(6,t) u(0,t)du(6,t) v d 2 u(6,t) 1 



-I{u(6,t)} 



(5) 



dt L 2 86 L 2 OO 2 

It is well known that the flat front solution of this equation is linearly unstable. The 
linear spectrum in fc-representation is 



oj k = \k\/L-vk 2 /L 2 



(6) 
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There exists a typical scale k max which is the last unstable mode 



L 



(7) 



Nonlinear effects stabilize a new steady-state which is discussed next. 

The outstanding feature of the solutions of this equation is the appearance of cusp- 
like structures in the developing fronts. Therefore a representation in terms of Fourier 
modes is very inefficient. Rather, it appears very worthwhile to represent such solu- 
tions in terms of sums of functions of poles in the complex plane. It will be shown 
below that the position of the cusp along the front is determined by the real coordinate 
of the pole, whereas the height of the cusp is in correspondence with the imaginary 
coordinate. Moreover, it will be seen that the dynamics of the developing front can 
be usefully described in terms of the dynamics of the poles. Following 1171 [T4l fj"5l [PD 
we expand the solutions u(9,t) in functions that depend on N poles whose position 
Zj(t) = Xj(t) + iyj(t) in the complex plane is time dependent: 



JV 

u(6, t) = v cot 



+ c.c. 



2 sin[0 - Xj {t)\ 



* cosh[j/j (t)] — cos[6> — Xj(t)] 



(8) 



JV 

h(6,t) = 2i/^ In [cosh^O)) - cos(0 - xj(t)) +C(t) 

3=1 



(9) 



In (0 C(t) is a function of time. The function (O is a superposition of quasi-cusps 
(i.e. cusps that are rounded at the tip). The real part of the pole position (i.e. xj) 
is the coordinate (in the domain [0, 2tt]) of the maximum of the quasi-cusp, and the 
imaginary part of the pole position (i.e yj) is related to the depth of the quasi-cusp. As 
yj decreases the depth of the cusp increases. As yj —> the depth diverges to infinity. 
Conversely, when yj — »• oo the depth decreases to zero. 

The main advantage of this representation is that the propagation and wrinkling of 
the front can be described via the dynamics of the poles. Substituting (|8]l in (|5]l we 
derive the following ordinary differential equations for the positions of the poles: 



L 



< dzj 
dt 



2JV 

■ E 



cot 



Zj - z k 



+ i—sign[Im(zj)} 



(10) 



We note that in (|8), due to the complex conjugation, we have 2A^ poles which are 
arranged in pairs such that for j < N Zj+n = %■ in the second sum in ((HJ each pair 
of poles contributed one term. In Eq.dTOb we again employ 2N poles since all of them 
interact. We can write the pole dynamics in terms of the real and imaginary parts Xj 
and yj. Because of the arrangement in pairs it is sufficient to write the equation for 
either yj > or for yj < 0. We opt for the first. The equations for the positions of the 
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poles read 



, dx 



N 



2 tlx j ^-—\ , 

L — = v s M x j - x k) 



[cosh(y j - y k ) 



(11) 



k=l,k^j 

■ cos(xj - Xk)}^ 1 + [cosh(y J + y k ) - cosixj - x k )]~ x 



dt 



N 

» E 

sinh(j/j + 



sinhfa- - y k ) 



vcosh^ 
Vk) 



Vk) 



- cos(xj — x k ) 
^coth(yj) — L 



(12) 



cosh(j/j + y k ) — 003(2;^ 

We note that if the initial conditions of the differential equation <(5j are expandable in 
a finite number of poles, these equations of motion preserve this number as a function 
of time. On the other hand, this may be an unstable situation for the partial differential 
equation, and noise can change the number of poles. This issue will be examined at 
length in Section [3] We turn now to a discussion of the steady state solution of the 
equations of the pole-dynamics. 



2.1 Qualitative properties of the stationary solution 

The steady-state solution of the flame front propagating in channels of width 2n was 
presented in Ref.[ 15 1. Using these results we can immediately translate the discussion 
to a channel of width L. The main results are summarized as follows: 

1. There is only one stable stationary solution which is geometrically represented 
by a giant cusp (or equivalently one finger) and analytically by N(L) poles which 
are aligned on one line parallel to the imaginary axis. The existence of this 
solution is made clearer with the following remarks. 

2. There exists an attraction between the poles along the real line. This is obvious 
from Eq.jTTb in which the sign of dxj/dt is always determined by sin(xj — x k ). 
The resulting dynamics merges all the x positions of poles whose y-position 
remains finite. 

3. The y positions are distinct, and the poles are aligned above each others in posi- 
tions yj~i < yj < yj+i with the maximal being yjv(L) ■ This can be understood 
from Eq.(fT2l in which the interaction is seen to be repulsive at short ranges, but 
changes sign at longer ranges. 

4. If one adds an additional pole to such a solution, this pole (or another) will be 
pushed to infinity along the imaginary axis. If the system has less than N(L) 
poles it is unstable to the addition of poles, and any noise will drive the system 
towards this unique state. The number N(L) is 



N(L) 



1 (L 
2U + 1 



(13) 
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where ... is the integer part. To see this consider a system with N poles and 
such that all the values of yj satisfy the condition < yj < y m ax- Add now one 
additional pole whose coordinates are z a = (x a ,y a ) with y a ^> y max - From the 
equation of motion for y a , (fT2l i we see that the terms in the sum are all of the 
order of unity as is also the cot(y a ) term. Thus the equation of motion of y a is 
approximately 

The fate of this pole depends on the number of other poles. If N is too large the 
pole will run to infinity, whereas if N is small the pole will be attracted towards 
the real axis. The condition for moving away to infinity is that N > N(L) where 
N(L) is given by (13[ . On the other hand the y coordinate of the poles cannot hit 
zero. Zero is a repulsive line, and poles are pushed away from zero with infinite 
velocity. To see this consider a pole whose yj approaches zero. For any finite 
L the term coth(i/j) grows unboundedly whereas all the other terms in Eq.dT2b 
remain bounded. 

5. The height of the cusp is proportional to L. The distribution of positions of the 
poles along the line of constant x was worked out in 1 1 5 1 . 

We will refer to the solution with all these properties as the Thual-Frisch-Henon (TFH)- 
cusp solution. 



2.2 Nonlinear Stability 

The intuition gained so far can be used to discuss the issue of stability of a stable system 
to larger perturbations. In other words, we may want to add to the system poles at finite 
values of y and ask about their fate. We first show in this subsection that poles whose 
initial y value is below y max ~ log(£ 2 /^ 2 ) will be attracted towards the real axis. The 
scenario is similar to the one described in the last paragraph. 

Suppose that we generate a stable system with a giant cusp at 8 C = with poles 
distributed along the y axis up to y max - We know that the sum of all the forces that 
act on the upper pole is zero. Consider then an additional pole inserted in the position 
(tt, y ma x)- It is obvious from Eq. dT2b that the forces acting on this pole will pull it 
downward. On the other hand if its initial position is much above y max the force on it 
will be repulsive towards infinity. We see that this simple argument identifies y max as 
the typical scale for nonlinear instability. 

Next we estimate y max and interpret our result in terms of the amplitude of a per- 
turbation of the flame front. We explained that uppermost pole's position fluctuates 
between a minimal value and infinity as L is changing. We want to estimate the char- 
acteristic scale of the minimal value of y max {L). To this aim we employ the result of 
ref.| 15 1 regarding the stable distribution of pole positions in a stable large system. The 
parametrization of lfl5l differs from ours; to go from our parametrization in Eq.© to 
theirs we need to rescale u by L^ 1 and t by L. The parameter v in their parameteri- 
zation is v/L in ours. According to [15] the number of poles between y and y + dy is 
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given by the p(y)dy where the density p(y) is 

p(y) = -^ln[coth(|y|/4)]. (15) 



To estimate the minimal value of y max we require that the tail of the distribution p(y) 
integrated between this value and infinity will allow one single pole. In other words, 

dyp(y) « 1 . (16) 

Expanding ( fT~5T > for large y and integrating explicitly the result in ([Tol l we end up with 
the estimate 

r 4L 

ymax ~ ^ In 



(17) 



For large L this result is y max ~ m (^2 )■ If we now add an additional pole in the 
position (9, y max ) this is equivalent to perturbing the solution u(9, t) with a function 
ve -Vmam s i n (0) 5 as can be seen directly from ([8). We thus conclude that the system is 
unstable to a perturbation larger than 

u{9) ~ v z sin(6>)/L 2 . (18) 

This indicates a very strong size dependence of the sensitivity of the giant cusp solution 
to external perturbations. This will be an important ingredient in our discussion of 
noisy systems. 



3 Acceleration of the Flame Front, Pole Dynamics and 
Noise 

A major motivation of this Section is the observation that in radial geometry the same 
equation of motion shows an acceleration of the flame front. The aim of this section is 
to argue that this phenomenon is caused by the noisy generation of new poles. More- 
over, it is our contention that a great deal can be learned about the acceleration in radial 
geometry by considering the effect of noise in channel growth. In Ref. lfl5l it was 
shown that any initial condition which is represented in poles goes to a unique station- 
ary state which is the giant cusp which propagates with a constant velocity v = 1/2 
up to small 1/L corrections. In light of our discussion of the last section we expect 
that any smooth enough initial condition will go to the same stationary state. Thus if 
there is no noise in the dynamics of a finite channel, no acceleration of the flame front 
is possible. What happens if we add noise to the system? 

For concreteness we introduce an additive white-noise term r}(9, t)lo the equation 
of motion (0 where 

77(^,0 =$> fc (i)exp(ifc0), (19) 

k 

and the Fourier amplitudes rjk are correlated according to 

< Vk(t)vl>{t') >= ts Kk ,5(t - t') . (20) 
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We will first examine the result of numerical simulations of noise-driven dynamics, and 
later return to the theoretical analysis. 

3.1 Noisy Simulations 

Previous numerical investigations j5] [23 1 did not introduce noise in a controlled fash- 
ion. We will argue later that some of the phenomena encountered in these simulations 
can be ascribed to the (uncontrolled) numerical noise. We performed numerical simu- 
lations of Eq. (fusing a pseudo-spectral method. The time-stepping scheme was chosen 
as Adams-Bashforth with 2nd order presicion in time. The additive white noise was 
generated in Fourier-space by choosing r/k for every k from a flat distribution in the 

interval [— \J 2j-, ij 2j-] . We examined the average steady state velocity of the front as 
a function of L for fixed / and as a function of / for fixed L. We found the interesting 
phenomena that are summarized here: 

1. In Fig Q] we can see two different regimes of behavior the average velocity v as 
function of noise / 5 for fixed system size L. For the noise / smaller then same 
fixed value f cr 

v~f*. (21) 

For these values of / this dependence is very weak, and £ w 0.02. For large 
values of / the dependence is much stronger 

2. In Fig|2]we can see growth of the average velocity v as function of the system 
size L. After some values of L we can see saturation of the velocity. For regime 
/ < fcr the growth of the velocity can be written as 

v~U\ /i w 0.40 ± 0.05 . (22) 

3.2 Calculation of Poles Number in the System 

The interesting problem that we would like to solve here it is to find number of poles 
that exist in our system outside the giant cusp. We can make it by next way: to calculate 
number of cusps (points of minimum or inflexional points) and their position on the 
interval 9 : [0, 2tt] in every moment of time and to draw positions of cusp like function 
of time, see Figj3] 

We assume that our system is almost all time in "quasi-stable" state, i.e. every new 
cusp that appears in the system includes only one pole. By help pictures obtained by 
such way we can find 

1 . By calculation number of cusp in some moment of time and by investigation of 
history of every cusp (except the giant cusp) , i.e. how many initial cusps take 
part in formation this cusp, after averaging with respect to different moments of 
time we can find mean number of poles that exist in our system outside the giant 
cusp. Let us denote this number SN. We can see four regimes that can be define 
with respect to dependence of this number on noise /: 
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(i) Regime I: Such small noise that no poles exist in our system outside the giant 
cusp. 

(ii) Regime II Strong dependence of poles number SN on noise /. 

(iii) Regime III Saturation poles number SN on noise /, so we see very small 
dependence of this number on noise 

SN ~ /°- 03 (23) 
Saturated value of SN is defined by next formula 

<5iV w N(L)/2 w i- (24) 
where N(L) w ~ is number of poles in giant cusp. 

(iv) Regime IV We again see strong dependence of poles number SN on noise 
/■ 

SN ~ / 01 (25) 

Because of numerical noise we can see in most of simulations only regime III 
and IV. In future if we don't note something different we discuss regime III. 



2. By calculation of new cusp number we can find number of poles that appear in 

dN 
dt 



the system in unit of time . In regime III 



dN #0.03 nf.\ 
-dt ~ 7 (26) 



Dependence on L and v define by 

dN 



L 0,8 (27) 
dt 



f ~ ^ < 28 » 



And in regime IV 

dN 



dt (29) 
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3.3 Theoretical Discussion of the Effect of Noise 



3.3.1 The Threshold of Instability to Added Noise. Transition from regime I to 
regime II 

First we present the theoretical arguments that explain the sensitivity of the giant cusp 
solution to the effect of added noise. This sensitivity increases dramatically with in- 
creasing the system size L. To see this we use again the relationship between the linear 
stability analysis and the pole dynamics. 

Our additive noise introduces perturbations with all /c-vectors. We showed previ- 
ously that the most unstable mode is the k = 1 component A\ sin(0). Thus the most 
effective noisy perturbation is rji sin(0) which can potentially lead to a growth of the 
most unstable mode. Whether or not this mode will grow depends on the amplitude of 
the noise. To see this clearly we return to the pole description. For small values of the 
amplitude A\ we represent A\ sin(0) as a single pole solution of the functional form 
i/e~ v sin 9. The y position is determined from y = — log |Al|/V, and the ^-position is 
9 = 7r for positive A\ and 9 — for negative A% . For very small A\ the fate of the pole 
is to be pushed to infinity, independently of its 9 position; the dynamics is symmetric in 
A\ — ^ — A\ when y is large enough. On the other hand when the value of A\ increases 
the symmetry is broken and the 9 position and the sign of A\ become very important. 
If A\ > there is a threshold value of y below which the pole is attracted down. On 
the other hand if A\ < 0, and 9 — the repulsion from the poles of the giant cusp 
grows with decreasing y. We thus understand that qualitatively speaking the dynamics 
of Ai is characterized by an asymmetric "potential" according to 

Al dXV ' (30) 

V(A{) = \A\-aA\ + ... . (31) 

From the linear stability analysis we know that A s» vj L 2 , cf. Eq.dTTb. We know 
further that the threshold for nonlinear instability is at A\ ps v z /L 2 , cf. EqdT8T>. This 
determines that value of the coefficient a pa 2/2>v 2 . The magnitude of the "potential" 
at the maximum is 

V(A max ) « z//L 6 . (32) 

The effect of the noise on the development of the mode Ai sin 9 can be understood 
from the following stochastic equation 

A, = -MM +,»«,. (33) 

It is well known [24| that for such dynamics the rate of escape R over the "potential" 
barrier for small noise is proportional to 

i?~^exp-" 7 /^ 6 . (34) 

The conclusion is that any arbitrarily tiny noise becomes effective when the system 
size increase and when v decreases. If we drive the system with noise of amplitude t 
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the system can always be sensitive to this noise when its size exceeds a critical value 
L c that is determined by -j- ~ v 7 /L®. This formula defines transition from regime 
I (no poles) to regime II. For L > L c the noise will introduce new poles into the 
system. Even numerical noise in simulations involving large size systems may have a 
macroscopic influence. 

The appearance of new poles must increase the velocity of the front. The velocity 
is proportional to the mean of (ujV) 2 . New poles distort the giant cusp by additional 
smaller cusps on the wings of the giant cusp, increasing u 2 . Upon increasing the noise 
amplitude more and more smaller cusps appear in the front, and inevitably the velocity 
increases. This phenomenon is discussed quantitatively in Section[3] 



3.3.2 Verifying of asymmetric "potential" form 

From the equations of the motion for poles we can find the distribution of poles in 
the giant cusp lfT5l . If we know the distribution of poles in the giant cusp we can 
then find the form of the "potential" and verify numerically expressions for values A, 
Amax and dV g^ discussed previously. The connection between amplitude A\ and 
the position of the pole y is defined by Ai — AveT v and the connection between 
the potential function dV QA^ an< ^ tne P ^ 011 °f tne P°l e V i s defined by formula 
= 4:V^e~ v , where ^ can be determined from the equation of the motion 
of the poles. We can find A max as the zero-point of 9 g^ 1 " 1 and A can be found as 

I 8 aA* 1 ^ ^ or ^ = 0' Numerical measurements were made for the set of values 
L = 2nv, where n is a integer and n > 2. For our numerical measurements we use 
the constant v = 0.005 and the variable L, where L changes in the interval [1,150], or 
variable v that changes in the interval [0.005,0.05] and the constant L = 1. The results 
obtained follow: 



1 . Formula for 



2. Formula for 4* 



An(l) 



6.5 . (35) 



0.465 . (36) 



where A N ^ is defined by position of the upper pole. 

\ T 2 

Formula for 

V 

XL 2 



= 0.5 . (37) 
v 



We also verify the boundary between regime I (no new cusps) and regime II 

J_ 



(new cusps appear). Fig. |4]shows the dependence of f- on L c . We can see that 



f/L c - 1/Ll 
These results are in good agreement with the theory. 
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3.3.3 The Noisy Steady State and its Collapse with Large Noise and System Size 

In this subsection we discuss the response of the giant cusp solution to noise levels that 
are able to introduce a large number of excess poles in addition to those existing in the 
giant cusp. We will denote the excess number of poles by SN. The first question that 
we address is how difficult is it to insert yet an additional pole when there is already 
a given excess SN. To this aim we estimate the effective potential Vsn{A\) which is 
similar to PTT ) but is taking into account the existence of an excess number of poles. 
A basic approximation that we employ is that the fundamental form of the giant cusp 
solution is not seriously modified by the existence of an excess number of poles. Of 
course this approximation breaks down quantitatively already with one excess pole. 
Qualitatively however it holds well until the excess number of poles is of the order of 
the original number N(L) of the giant cusp solution. Another approximation is that the 
rest of the linear modes play no role in this case. At this point we limit the discussion 
therefore to the situation SN <C N(L) (regime II). 

To estimate the parameter A in the effective potential we consider the dynamics of 
one pole whose y position y a is far above y m ax- According to Eq.(fT4b the dynamics 

rCadS dy a ^ 2u(N(L) + SN) 1 

~dT T? I (38) 

Since the N(L) term cancels against the L _1 term (cf. Sec. 12. U . we remain with a 
repulsive term that in the effective potential translates to 

vSN 

X = — - ( 39 > 

Next we estimate the value of the potential at the break-even point between attraction 
and repulsion. In the last subsection we saw that a foreign pole has to be inserted below 
ymax in order to be attracted towards the real axis. Now we need to push the new pole 
below the position of the existing pole whose index is N(L) — SN. This position is 
estimated as in Sec l2.2l bv employing the TFH distribution function (1151 1. We find 



ysN ~ 2 In 



4L 



;5N 



(40) 



As before, this implies a threshold value of the amplitude of single pole solution 
Amax sin 6 which is obtained from equating A max — ve~ V5N . We thus find in the 
present case A max ~ v i (SN) 2 / L 2 . Using again a cubic representation for the effec- 
tive potential we find a = 2/ {3v 2 SN) and 

1 v 7 (SN) 5 

V(A max ) = g^r^ ■ (41) 

Repeating the calculation of the escape rate over the potential barrier we find in the 
present case 

i^^exp-W/^ 6 . (42) 

For a given noise amplitude t there is always a value of L and v for which the 
escape rate is of O(l) as long as SN is not too large. When SN increases the escape 
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rate decreases, and eventually no additional poles can creep into the system. The typical 
number SN for fixed values of the parameters is estimated from equating the argument 
in the exponent to unity: 

M \ 1/5 
SN w f jL*fv 7 \ . (43) 

We can see that SN depend on noise / very seriously. It is not the case in regime 
III. Let us find conditions of transition from regime II to regime III, where we see 
saturation of SN with respect to noise /. 

(i) We use for the amplitude of pole solution that really equal to cosh (y s s ^l cosg 
expression A max = Ave~ V5N but it is right only for big number ysN- For y$N < 

1 better approximation is A rnax = -4r-. From equation (l40l > we can find that the 

y&N 

boundary value ysN = 1 correspond to SN s» N(L)/2 



(ii) We use expression ysN ~ 2 In 



4L 



but for big value of SN better approxi- 



8eL 

tt' 2 v{N(L)-8N) 



mation that can be find the same way is ygN ~ ^T" (N(L) — SN) In 
These expressions give us nearly equal result for SN ps N(L)/2. 
From (i) and (ii) we can make next conclusions 

(a) Transition from regime II to regime III happens for nearly SN ~ N(L) /2 

(b) Using new expression in (i) and (ii) for amplitude A max and ysN we can find 
for noise in regime III: 

L - V[A max ) ~ \A max ~ ^ (^) ~ ^ ( N{L) _ SN) 4 ( 44 > 

This expression define very slow dependence of SN on noise j- for SN > N(L) /2 
that explain noise saturation of SN for regime III. 

(c) Form of the giant cusp solution is defined by poles that are closely to zero with 
respect toy. Forregime III N(L)/2 poles that have position y < ysN=N(L)/2 = 1 sta y 
on these place. This result explain why giant cusp solution is not seriously modified 
for regime III. 

From Eq.d43b by help of boundary condition 

SN w N(L)/2 (45) 
boundary noise between regime II and III can be found 

h - v 2 (46) 
The basic equation describing pole dynamics is next 

where ^ is number of poles that appear in unit of time in our system, SN is excess 
number of poles, T is mean life time of pole (between appearing and merging with giant 
cusp). Using result of numerical simulations for ^ and d45t we can find for T 
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T = § ~ ^ (48) 

dt 

So life time proportional to v and depend on system size L very weakly. 

Moreover, the lifetime of a pole is defined by the lifetime of the poles that are in 
a cusp. From the maximum point of the linear part of Eq.(Q]), we can find the mean 
character size 

X m ~ v (49) 

that defines the size of our cusps. This result was confirmed in numerical calcula- 
tions executed in 0251 , Indeed, we can see it from Fig. 9 in [25]. The mean number of 
poles in a cusp 

n big « ~ const (50) 
does not depend on L and v. The mean number of cusps is 

N big ~ ™- ~ L - . (51) 

nu g v 

Let us assume that some cusp exists in the main minimum of the system. The 
lifetime of a pole in such a cusp is defined by three parts. 

(I) Time of the cusp formation. This time is proportional to the cusp size (with 
ln-corrections) and the pole number in the cusp (from pole motion equations) 

Ti ~ X m n b ig ~ v (52) 

(II) Time that the cusp is in the minimum neighborhood. This time is defined by 

T 2 ~ - (53) 

v 

where a is a neighborhood of minimum, such that the force from the giant cusp is 
smaller than the force from the fluctuations of the excess pole number SN, and v is the 
velocity of a pole in this neighborhood. Fluctuations of excess pole number SN are 
expressed as 

N fl = VSN . (54) 
From this result and the pole motion equations we find that 

The velocity from the giant cusp is defined by 
So from equating these two equations we obtain 
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Thus for T2 we obtain 



T 2 ~ - ~ L . (58) 

v 

(III) Time of attraction to the giant cusp. From the equations of motion for the 
poles we get 

T 3 ~ Lln(-) - LlnVX~ L . (59) 
a 

The investigated domain of the system size was found to be 

Ti»T 2 ,T 3 (60) 

Therefore full lifetime is 

T = T 1 +T 2 +T 3 ~v + S L, (61) 
where s is a constant and 

< s < 1 . (62) 

This result qualitatively and partly quantatively explains dependence d48l . From 
d48l . ( f47b , (l45l l we can see that in regime III ^ is saturated with the system size L. 

3.4 The acceleration of the flame front due to noise 

In this section we estimate the scaling exponents that characterize the velocity of the 
flame front as a function of the system size. Our arguments in this section are even less 
solid than the previous ones, but nevertheless we believe that we succeed to capture 
some of the essential qualitative physics that underlies the interaction between noise 
and instability and which results in the acceleration of the flame front. 

To estimate the velocity of the flame front we need to write down an equation for 
the mean of < dh/dt > given an arbitrary number TV of poles in the system. This 
equation follows directly from ©: 



dh\ 1 1 ^ 



at I L z 2n J Q 
After substitution of © in d63l we get, using (TTTb and (fT2b 

1 fc=i v 7 

The estimates of the second and third terms in this equation are straightforward. Writ- 
ing N = N(L) + SN(L) and remembering that N(L) ~ L/v and SN(L) - N(L)/2 
we find that these terms contribute 0(1). The first term will contribute only when 
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the current of poles is asymmetric. Since noise introduces poles at a finite value of 
Vmin, whereas the rejected poles stream towards infinity and disappear at boundary of 
nonlinearity defined by position of highest pole 

r 4£ i 

Vmax ~ 2 In . (65) 

, we have an asymmetry that contributes to the velocity of the front. To estimate 
the first term let us define 

, l+dl , 

where iMK [ s sum over poles that are on the interval y : [I, I + dl]. We can 

write 



ing down. 

For these flux we can write 



dt ' y ^ dt ' ap dt 
Where d(J2 -J^)«p A ux of poles moving up and d(J2 ^ d j L )down flux of poles mov- 



d{}^-7T)up,d(2^-7-)down< -7-dl (68) 



dt dt > auwn - dt 

So for the first term 



N A 



dl 



dl 



L 



dl 



dN 

~dt ax y m ^ n y 

dN_ 

— TrV max 

dt 

Because of slow(ln) dependence of y max on L and v ^ term define oder of non- 
linearity for first term. This term zero for symmetric current of poles and achieves 
maximum for maximal asymmetric current of poles. Comparison v ~ L aA2 /°' 02 and 

^ ^°- 8 / - 03 confirm this calculation. 
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4 summary and conclusions 



The main two messages of this paper like the previous one are: (i) There is an im- 
portant interaction between the instability of developing fronts and random noise; (ii) 
This interaction and its implications can be understood qualitatively and sometimes 
quantitatively using the description in terms of complex poles. 

The pole description is natural in this context firstly because it provides an exact 
(and effective) representation of the steady state without noise. Once one succeeds to 
describe also the perturbations about this steady state in terms of poles, one achieves 
a particularly transparent language for the study of the interplay between noise and 
instability. This language also allows us to describe in qualitative and semi-quantitative 
terms the inverse cascade process of increasing typical lengths when the system relaxes 
to the steady state from small, random initial conditions. 

The main conceptual steps in this paper are as follows: firstly one realizes that the 
steady state solution, which is characterized by N(L) poles aligned along the imag- 
inary axis is marginally stable against noise in a periodic array of L values. For all 
values of L the steady state is nonlinearly unstable against noise. The main and fore- 
most effect of noise of a given amplitude / is to introduce an excess number of poles 
5N(L, /) into the system. The existence of this excess number of poles is responsible 
for the additional wrinkling of the flame front on top of the giant cusp, and for the 
observed acceleration of the flame front. By considering the noisy appearance of new 
poles we rationalize the observed scaling laws as a function of the noise amplitude and 
the system size. 

The "phase diagram" as a function of L and / in this system consists of four 
regimes (in contradiction with our previous results |6|). In the first one, discussed 
in Section [3.3.1l . the noise is too small to have any effect on the giant cusp solution. 
The second regime (very small excess number of poles ) can not be observed because 
of numerical noise and discussed only theoretically. In the third regime the noise intro- 
duces excess poles that serve to decorate the giant cusp with side cusps. In this regime 
we find scaling laws for the velocity as a function of L and / and we are reasonably 
successful in understanding the scaling exponents. In the fourth regime the noise is 
large enough to create small scale structures that are not neatly understood in terms of 
individual poles. It appears from our numerics that in this regime the roughening of 
the flame front gains a contribution from the the small scale structure in a way that is 
reminiscent of stable, noise driven growth models like the Kardar-Parisi-Zhang model. 

One of our main motivations in this research was to understand the phenomena ob- 
served in radial geometry with expanding flame fronts. A full analysis of this problem 
cannot be presented here. We note however that many of the insights offered above 
translate immediately to that problem. Indeed, in radial geometry the flame front ac- 
celerates and cusps multiply and form a hierarchic structure as time progresses. Since 
the radius (and the typical scale) increase in this system all the time, new poles will be 
added to the system even by a vanishingly small noise. The marginal stability found 
above holds also in this case, and the system will allow the introduction of excess poles 
as a result of noise. The results discussed in Ref.[7| can be combined with the present 
insights to provide a theory of radial growth. This theory was offered in Ref.JS). 

We have had a serious open problem for this case [10], but this problem was solved 
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successfully recently Karlin and Sivashinsky 1 12| , lfl3l . For a cylindrical case of the 
flame front propagation problem at absence of noise (only numerical noise Ref. li26l 
27, 28 29 30]) by Sivashinsky with help of numerical methods it was shown, that the 
flame front is continuously accelerated. During all this account time it is not visible 
any attributes of saturation. To increase time of the account is a difficult task. Hence, 
absence or presence of velocity saturation in a cylindrical case, as consequence of the 
flame front motion equation it was a open problem before appearance of the papers 

02, 021 

For the best understanding of dependence of flame front velocity as functions of its 
radius in a cylindrical case similar dependence of flame front velocity on width of the 
channel (in a flat case) also was analyzed by numerical methods. Growth of velocity is 
also observed and at absence of noise (only numerical noise!) also any saturation of the 
velocity it is not observed. Introduction obvious Gaussian noise results to appearance 
of a point of saturation and its removal from the origin of coordinates with decreasing 
of noise amplitude, allowing extrapolating results on small numerical noise. (Figff]) 

Hence, introducing of Gaussian noise in numerical calculation also for a cylindrical 
case can again results to appearance of a saturation point and will allow to investigate 
its behavior as function of noise amplitude by extrapolating results on small numerical 
noise. This investigation was really made and saturation was observed recently by 
Karlin and Sivashinsky lfl2l , iTHl for ID, 2D and 3D cases. 

Finally, the success of this approach in the case of flame propagation demonstrates 
that Laplacian growth patterns can be dealt with using similar ideas. A problem of im- 
mediate interest is Laplacian growth in channels, in which a finger steady-state solution 
is known to exist. It is documented that the stability of such a finger solution to noise 
decreases rapidly with increasing the channel width. In addition, it is understood that 
noise brings about additional geometric features on top of the finger. There are enough 
similarities here to indicate that a careful analysis of the analytic theory may shed as 
much light on that problem as on the present one. 
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Figures Legends 

Fig.l: The dependence of the average velocity v on the noise / 5 for L=10, 40, 80. 

Fig. 2: The dependence of the average velocity v on the system size L for /° 5 = 
0, 2.7 * 10~ 6 , 2.7 * 10~ 5 , 2.7 * 10~ 4 , 2.7 * 10~ 3 , 2.7 * 10~ 2 , 2.7 * 10"\ 0.5, 1.3, 2.7 . 

Fig. 3: The dependence of the cusps positions on time.L = 80 v = 0.1 / = 9 * 10—6 

Fig. 4: The first odd eigenfunction obtained from traditional stability analysis. 
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